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Abstract 

We have studied the reaction rjd — > r/d close to threshold within a nonrela- 
tivistic three-body formalism. We considered several rjN and NN models, in 
particular potentials with separable form, fitted to the low-energy r/N and 
NN data to represent the two-body interactions. We found that with real- 
istic two-body interactions a quasibound state does not exist in this system, 
although there is an enhancement of the cross section by one order of mag- 
nitude, in the region near threshold, which is a genuine three-body effect not 

predicted within the impulse approximation. 
14.40.Ag,25.40.Ve,25.80.Hp 
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I. INTRODUCTION 



The elastic scattering r]d reaction has been studied recently by several authors in 
order to investigate the existence of a resonance or a quasibound state in this system, for 
which there are experimental indications @-^|. Some of these studies concluded that such 
an state would exist for certain values of the two-body rjN data. However, since they used in 
one form or another incomplete information on the two-body subsystems (in particular that 
corresponding to the rjN sector), we believe that a new calculation is required which takes 
into account all the information that is now available. In particular, we study the effect of 
the repulsion at short distances of the NN interaction and take into account the rjN — > rjN 
scattering amplitude that has been determined recently P-|T!|. 

We will present the Faddeev formalism in section II. In section III we will first calculate 
the eta-deuteron scattering length to compare with the results of multiple-scattering theories 
as well as with separable potentials models and finally we will present the predictions of our 
model for rjd scattering. We will give our conclusions in section IV. 

II. FADDEEV FORMALISM 

Let us consider a system of three particles, where two of them are identical, interacting 
pairwise through separable potentials that act only in S-waves. In the case of the r]d system, 
S-wave means, for the eta-nucleon pair the Sn channel, and for the nucleon-nucleon pair the 
3 Si channel. The two-body T-matrix of the pair jk will be assumed of the separable form 

kiPhPi, E ) = 9i(Pi)Ti(E)9M, (!) 

where Ti(E) and gi(pi) will be specified later. In the following we will identify particle 1 
with the rj and the identical particles 2 and 3 with the two nucleons. 

The Faddeev equations for the case of rjd scattering can be solved explicitely for the 
7]d — > N(Nrj) transition amplitude T2, with nucleon 2 being the spectator particle in the 
final state. One obtains 
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T 2 (q 2 ; E) = 2K 21 (q 2 , q w ; E) + / q' 2 2 dq' 2 M(q 2 , q' 2 ; E)r 2 (E - q' 2 /2v 2 )T 2 (q' 2] E), (2) 



o 



where 



M(q 2 ,q 2 ;E) = K 23 (q 2 ,q' 2 ;E) + 2 / q 2 1 dq 1 K 21 (q 2 ,q 1 ; E)n(E - q 2 1 /2u 1 )K 12 (q 1 ,q' 2 - E), (3) 

Jo 

and the driving term for the transition from the amplitude with spectator j in the initial 
state to the one with spectator i in the final state is 

1 f 1 A^a 9i{Pi)9j{Pj) 



with 



Pi = ( — 2 £ + + 2^q iqj cose] , (5) 
\mi J m k J 



/ 2 \ 1/2 

Pj = + 9< + ^—mjcose . (6) 



m% - m k 
fj,i and Vi are the reduced masses 



Hi = ~ , (7) 

rrij + m k 



mi(mj + m k ) 

^ = , 7 , 8 

m, + rrij + m k 

and the rjd on-shell relative momentum q w is defined by the relation 

where Bd is the binding energy of the deuteron. 

The r]d elastic-scattering amplitude is obtained from the solution of Eq. (|2]) by inserting 
into it a final rjd state, as 

jrjy-, /*oo 

F vd {qio) = — j^- J Q q 2 2 dq 2 K 12 (q 10 , q 2 ; E)r 2 (E - q 2 2 /2u 2 )T 2 (q 2 ; E). (10) 
where N is the normalization of the deuteron wave function 



3 



N — / v\dVi 



B d + p\ /2/ii 



(11) 



The r/d scattering length is given by 



(12) 



while the integrated elastic cross section is given by 



oelas = 47r|F^(^io) | 2 . 



(13) 



Notice that Eqs. (Q) - @ and fllQl) do not include the irN channel explicitely, but only 
through the inelasticity of the r)N channel. As for the the i]N inelasticity due to the nuN 
channel, its contribution is not yet included at this stage of the calculations. 



III. RESULTS 

We started by calculating the r/d scattering length (0), by solving the integral equation 
(0) with the method of matrix inversion after replacing the integration with a gaussian 
mesh. Notice that in this case Eqs. (^)-(^) are free of singularities ( the singularity of T\ at 
E = —Bd in eq. (|3]) occurs for qi = and is regularized by the integration volume element). 
Additionally, we also calculated the integrated elastic cross section of rjd scattering. In order 
to solve the integral equation (H) above threshold, we used the method of contour rotation 

m 



A. A^d with non- dynamical separable models 

We will calculate here the r]d scattering length A n ^ for the models proposed in Refs. 
The signal that a quasibound state exists for a given model is that the real part of 
A v d becomes negative while the imaginary part gets large. 

Notice that in Eq. (|I]) we have assumed a separable model for the two-body amplitudes ti. 
This form of the T-matrix is obtained if one assumes a separable potential between particles 
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j and k: a dynamical two-body equation determines the function Ti(E) by the form factors 
gi(pi), which carry information on the the range of the potential, and the strength parameter 
of the potential. However, in the AGS formalism |14| used in Ref. and in the multiple- 



scattering approach used in Ref. the function r%{E) for the 7]N subsystem, instead of being 
calculated, has been chosen independently of the form factor g^ip-i)- Such an assumption 
violates the spirit of the Faddeev approach which requires a two-body interaction in order to 
relate the off-shell behavior of the T-matrix in the energy variable E to the off-shell behavior 
in the momentum variable p2- Nevertheless, it is instructive to repeat those calculations in 
order to check the accuracy of our numerical solution by comparing with the exact result of 
Ref. as well as to test the convergence of the multiple-scattering schemme developed in 
which is based in a partial summation of the multiple-scattering series. 

In both Refs. the form factor (72(^2) has been taken of the Yamaguchi form 

= ^2—^2' ( 14 ) 

with « 2 = 3.316 fm~l. In Ref. ||, the function t 2 (E) has been parametrized as 

«w = E^kw (15) 

with Eq = 1535 MeV - (mjv + Tn v ) and T = 150 MeV. The parameter A,, in Eq. (|15D was 
chosen to reproduce the complex r/N scattering length a V N, by using the relation 

t 2 (0,0;0) = -^, (16) 

where \x-i is the rjN reduced mass. This leads to 

As for the nucleon-nucleon separable T-matrix used in Ref. it was generated from a 
Yamaguchi separable potential with an energy-dependent strength [||,[n|]. Using these pa- 



rameters we calculated the 7]d scattering length A v d with the formalism described in the 
previous section for a variety of values of the r]N scattering length a V N that have been pro- 
posed in the literature (see Refs. ||[4j]). We show the results of this comparison in table 



I. The results of Ref. [|J using the AGS formalism are given in column two and the ones 
obtained with the formalism of the previous section are given in column three. As one can 
see from table I there is very good agreement between the two calculations. The small dis- 
crepancies shown in some cases are of no significance since they occur for the values of a V N 
allowing the quasibound state to occur and thereby the solutions are highly unstable. In 
our case, in this situation we had to use a large number of mesh points in order to guarantee 
stability. 

In the multiple-scattering approach of Ref. || an approximate formula was used which is 
based in a partial sumation of the multiple-scattering series. The function T2(E) was taken 
to be constant 

t 2 (E) = X v . (18) 

Using the relation (|16|) this gives 

A„ = —a V N, (19) 

which will be referred to as their model I. They used also a second model which will be 
referred to as model II for which instead of Eq. ([TJ]) they took 

A, = (20) 

with go = i0.367 fin -1 . For the nucleon-nucleon interaction they used a Yamaguchi separable 
potential with a range parameter a\ = 1.41 fm~l. 

We compare in table II the results of our exact calculations which we obtained using the 
parameters of Ref. || with their results using an approximate formula for the two models I 
and II. As it can be seen from this table, the approximate formula of the multiple scattering 
series works very well for small values of expected from convergence arguments. 

When a V N is large the multiple scattering series formula deviates more from the exact result, 
nevertheless, it is still qualitatively correct, since it predicts correctly the quasibound states 
in all the cases where they exist for both models. 
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B. A^d with separable-potential models 

In the previous subsection we have seen that the models of Refs. predict a quasi- 
bound state if the real part of a V N is of the order of 0.7-0.8 fm. However, since their rjN 
T- matrix is not derived from a potential their function T2(E) is not constrained by their 
form factor g 2 {jp 2 )- We will therefore construct separable potential models of the coupled 
rjN - tcN system that reproduce in one case just the complex scattering length a V N for ar- 
bitrary values of the rjN range-parameter a 2 and in another case the full 7]N-i]N scattering 
amplitude around the Sn resonance. 

Similarly, we will consider two different models of the NN interaction; a simple Yam- 
aguchi potential that does not have short-range repulsion and a PEST model which has the 
same half-off-shell behavior as the Paris potential so that it contains short-range repulsion. 

If we use in the Lippmann-Schwinger equation of the coupled r/N - ttN system the 
separable model 

< p\V m \p' >= \92(p)g2(p'), (21) 



< p\Vitk\p' >= KgAp)9n{p') 



< p\V m \p' >= ±JX v X 7r g 2 (p)g n (p' 



the T-matrices are of the form 



<p\t m (E)\p' >= g 2 (p)r 2 (E)g 2 (p') 



< p\t ww (E)\p' >= ^ L g 7r (p)r 2 (E)g 7r (p) 



(22) 
(23) 



(24) 
(25) 



< p\t m (E)\p >= ± 



\ 



(26) 



with 
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CO 



G2(E) = I v 2 d P = - . , (28) 

Jo E — p 2 /2fi 2 + le 



GJE) = I p 2 dp- . (29) 

" V ; Jo y y E + p 2 /2^-p 2 /2^ + te V ; 

fi2 and /j, n are the ^iV and ttN reduced masses respectively while p is the ttN relative 
momentum at the r)N threshold, i.e., 

2 _ [go ~ (^tt + ^v) 2 ][gp ~ (m n - m N ) 2 ] 
Po- , W 

with 

s = (m, + m^) 2 . (31) 
If we use simple Yamaguchi form factors 

92(P) = (32) 

^tt(p) = 2 - 2 , (33) 

we find that the strengths X v and A„- can be obtained in terms of the real and imaginary 
parts of a V N as 

! _ _ a l - Po ^/fg Im a vN _ 7r/i 2 -Re 

v 2a\ 2a w p a\ \a vN \ 2 a\ \a vN \ 2 ' 



x V>2 (al+pl) 2 Ima vN 

Att — 4 j [7T- I -30 J 



Since we do not include the pion channel explicitly but only through the function t 2 (E) (see 
Eq. (p?D), we will fix the range of the nN potential to the value a n = po, for which case the 
second term in the r.h.s. of Eq. drops out and it becomes clear that in this case the 
strength of the rjN potential is determined by the real part of the rjN scattering length and 
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the strength of the ttN potential is determined by the imaginary part of the rjN scattering 
length for a given value of the range a%. Since this models are based in a Yamaguchi form 
factor for the rjN potential we will refer to them as Y v n models. 

We constructed also more realistic separable models that reproduce not only the rjN 
scattering length but also the most important features of the Sn resonance such as its 
position and width. For this we considered the Sn amplitudes obtained from the analyses 
of Refs. [[ToH l^l. We found that with a simple Yamaguchi model of the rjN form factor is 
not possible to generate a resonance in the i] N Sn channel. We therefore changed the form 
factor (?2(p) instead of Eq. (^) to 

M = (^T7p> ( 36 ) 



while keeping for the tcN form factor the Yamaguchi form (|33|). We give in table III the 
parameters a n , \ n , «2, A, X v of the coupled rjN-irN separable potentials fitted to the Sn 
amplitudes of [p70|JTT|| as well as to the models A, B, C, and D of ||12|| . We show in Fig. 1, 



as an example, the r/N-r/N amplitude of Ref. |L0| (dashed lines) compared with the ones 
of our separable-potential model. Similar results are obtained for the other models. Since 
these models generate a resonance in the r]N channel we will refer to them as R^m models. 
Since we do not include the pion channel explicitly, only the t m component of the coupled 
i]N-ttN T-matrix given by Eq. (^) has been used after identifying t m with £2 of Eq. (1). 

In the case of the NN interaction we have considered two models; the simple Yamaguchi 
model used in Ref. || which has a range parameter ot\ = 1.41 fm _1 (which we will refer to 
as the Ynn model) and the PEST potential constructed in Ref. |K| (which we will refer to 
as the Pnn model) that is of the form 

6 £j 

0i(Pi) = E 3tV' ( 3? ) 

n=l Pi + P n 



where the parameters C n and (3 n are given in Ref. ||16| . The half-off-shell T-matrix of this 



separable potential has the same behavior as that of the Paris potential and therefore it 
takes into account the repulsion at short distances that is present in the nucleon-nucleon 
force. 

9 



We give in table IV the results of our separable-potential models for the r\d scattering 
length A V d where we have considered all four combinations of the rjN and NN separable- 
potential models. In the case of the rjN Yamaguchi model Y v n we took the range parameter 
a?2 = 3.316 fm -1 which is the same as in Refs. @||. Here however the corresponding T- 
matrix is calculated through the Lippmann-Schwinger equation. In the first column we 
give the reference for the rjN Sn amplitude that we used to fit the R v n model and in the 
second column we give the rjN scattering length of that amplitude which has been used to 
construct the Y v n model. The third column gives A v d using simple Yamaguchi models for 
the rjN and NN interactions and it shows that even with these simple separable models 
the quasibound state only appears when Rea^ is about 1.05-1.07 fm while in the multiple- 
scattering approaches of the previous section it appeared already with Rea V N ~ 0.6-0.9 fm. 
The fourth column gives the results of the Yamaguchi model for the rjN amplitude and the 
PEST model for the NN amplitude and it shows that the NN short-range repulsion also 
works against quasibinding since it wipes out the quasibound state although the ImA r]d 
remains large. The fifth and sixth columns contain the results of the resonant model of the 
rjN amplitude with Yamaguchi and PEST models for the NN interaction respectively, and 
they show that the attraction of the system is greatly reduced when one takes into account 
the resonant nature of the rjN amplitude. Notice, however, that experimentally the effects 
will be similar whether there is a quasibound state or not since in both cases there will be 
an enhancement of the cross section at threshold. 

C. rjd scattering 

We calculated the integrated elastic cross section of rjd scattering in the region of the 
Sn resonance for the six resonant models of the rjN interaction given in table III and the 
realistic PEST potential for the NN interaction. We show in Fig. 2 the results of the three- 
body model (solid lines) and of the impulse approximation (dashed lines). At threshold, 
the results of all the three-body models are about one order of magnitude larger than those 
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of the impulse approximation while at higher energies they are 2 or 3 times smaller. The 
behavior of the cross section at threshold indicates that even though the quasibound state 
is not present, the interaction in this region is very strong since it enhances the impulse 
approximation result by about one order of magnitude. Thus, a signal of this behavior may 
appear also in other processes where there is an nN final state like the np — > nd reaction 
where a large enhancement in the cross section has been observed in the region near threshold 

0- 

In order to illustrate the effect of the strong nd interaction in the reaction np — > nd 
we will estimate the enhancement of the np — > nd cross section at threshold due to the nd 
rescattering. We write the amplitude of the process np — > nd as 

A = B + BG T vd7 (38) 

where B is the amplitude of the np — > nd process without nd rescattering, Go is the two-body 
Lippmann-Schwinger propagator of the intermediate nd state, and T vd is the half-off-shell 
T-matrix of the elastic rjd process. If we introduce the nd elastic-scattering amplitude 
F vd = —7Tp: vd T vd , where fi vd is the rjd reduced mass, then at threshold, Eq. fl3"8] ) is written 
explicitly as 

A = B(0)[l + -J o d qi ^F vd ( gi )}, (39) 



where F v d(qi) is given by Eq. (|T0|) with qi replaced by q\. Therefore, the enhancement 



factor of the np — > nd cross section due to nd rescattering is 

2 r.,, B Mv ,.. m2 



/=|1+ w„ d *im FM ^ m 

The amplitude B of the np — > nd process without nd rescattering is presumably given 
by meson exchanges such as it, p, and n followed by the excitation and decay of the Sn 
resonance |[17|| . The explicit form of the production operator without nd rescattering is not 
so important since in our estimate of the enhancement factor given by Eq. ( PH) only the 
ratio B(qi)/B(0) enters. Therefore, we take for it the 77-exchange amplitude generated by 
our three-body model, i.e., 
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B(qi) = / q 2 2 dq 2 K l2 (qi, q 2 ; E)r 2 {E - q 2 2 /2v 2 )K 23 (q 2 , q 3 ; E), (41) 
Jo 

where E = —Bd and q% = yj (m v + mi) 2 — m 2 N is the momentum corresponding to an initial 
NN state. Using the models 1-6 of Fig. 2, taken from references [|T0|-[T2l and included in 
Table IV, we obtained the enhancement factors / = 2.5, 2.7, 3.1, 3.3, 4.7, and 5.1 respectively. 
These estimates are quite comparable to the enhancement factors observed in Ref. 0, in 
special for the 4 first models. 



IV. CONCLUSIONS 



Recently, the experimental band for the rjN scattering length a V N has been pushed 



towards larger values for its real part |Tl|Jr2|] . For those larger values, the rjN models not 
generated directly from an integral equation, which would fix naturally their off-energy- 
shell behavior needed in three-body calculations, predict a quasi-bound r]NN state. If 
instead, the new data is used to generate wN t-matrices calculated from a potential and 
an integral equation, our results indicate that a realistic NN interaction, like the Paris 
potential, through its short-range correlations, prevents the existence of the bound-state, 
independently of the r/N models, provided they have been built dynamically. We confirmed 
then that the predictions of a wNN eta-mesic nucleus are crucially affected by the off-shell 
behavior of the underlying n-N models. Importantly, however, is that even for an inexisting 
quasi-bound state, an exact three-body calculation for the multi-scattering series in the final 
state predicts a severe enhancement of the elastic rjd cross-section in the narrow region from 
threshold to 5-10 MeV above threshold. This result is independent of the 7]N two-body 
models. Very likely, the enhancement predicted by the exact three body calculations is 
related to the one observed in the reactions np — >• rjd and jd — > rjd HQ. We actually 
made an estimate of the enhancement of the cross section of np —>■ rjd due to the rjd final 
state interaction. Within the three-body model of our work, this enhancement is in the ball 
park of the empirical findings of 0] . 
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FIGURES 

FIG. 1. Real and imaginary parts of the r]N — ► r]N Sn amplitude. The solid and dashed lines 
are the results of our resonant R^n separable potential model 1; the symbols give the amplitudes 
ofRef. [10]. 

FIG. 2. Integrated elastic cross sections of the three-body model (solid lines) and of the im- 
pulse approximation (dashed lines) for the six resonant r/N models of table III using for the NN 
interaction the PEST model, as a function of the cm. rjd kinetic energy. 
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TABLES 

TABLE I. r/d scattering length (in fm) as calculated in Ref. [4] with the AGS formalism 
and as calculated here with the Faddeev formalism for different values of the rjN scattering length 
a v N (in fm). 



AGS [4] 



Faddeev 



0.25+i0.16 
0.30+10.30 
0.46+10.29 
0.579+i0.399 
0.876+10.274 
0.98+i0.37 



0.73+i0.56 
0.61+il.22 
1.31+11.99 
0.34+i3.31 
-8.81+14.30 
-4.69+11.59 



0.73+i0.56 
0.61+U.22 
1.33+11.99 
0.34+i3.33 
-9.02+14.15 
-4.74+11.53 



TABLE II. rjd scattering length (in fm) for models I and II calculated in Ref. [3] using a 
multiple-scattering theory (MST) and as calculated here using the Faddeev formalism for different 
values of the rjN scattering length a v N (in fm) . 



a V N 


MST I [3] 


Faddeev 


MST II [3] 


Faddeev 


0.25+i0.16 


0.66+i0.71 


0.64+i0.71 


0.66+i0.58 


0.65+i0.58 


0.30+i0.30 


0.39+il.28 


0.38+il.25 


0.58+il.ll 


0.56+H.09 


0.46+10.29 


0.72+i2.04 


0.62+11.95 


l.ll+il.54 


1.04+H.54 


0.579+i0.399 


-0.13+i2.64 


-0.08+i2.31 


0.93+i2.41 


0.74+i2.28 


0.876+i0.274 


-2.76+i4.24 


-1.54+i2.55 


2.42+15.55 


0.46+i5.21 


0.98+10.37 


-2.75+12.77 


-1.16+12.05 


-0.06+16.20 


-1.12+i4.21 
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TABLE III. Parameters of the t]N-ttN separable potential models fitted to the Sn resonant 
amplitudes given in Refs. [10-12]. 



Model 


Ref. 


a V N 




X n 


C*2 


A 


A^ 


1 


[6] 


0.72+10.26 


1.2 


-0.028013 


15.5 


21.394872 


-8817.102932 


2 


[7] 


0.75 + 10.27 


1.2 


-0.028013 


14.5 


18.657140 


-7222.670614 


3 


[8](D) 


0.83 + i0.27 


1.2 


-0.028013 


13.68 


17.7754 


-5970.573595 


4 


[8] (A) 


0.87+10.27 


1.2 


-0.028013 


13.1 


16.646837 


-5215.563508 


5 


[8](B) 


1.05 + 10.27 


1.2 


-0.028013 


11.1 


13.357879 


-3080.397824 


6 


[8](C) 


1.07 + i0.26 


1.2 


-0.028013 


10.47 


11.952401 


-2579.494411 



TABLE IV. rjd scattering length (in fm) for various models of the i]N and NN interactions 
where the respective t-matrices were calculated through a dynamical equation. Y v n stands for the 
Yamaguchi model of the rjN interaction with a range parameter 012 = 3.316 fm -1 , R v n stands 
for the resonant model of the t]N interaction, Ynn stands for the Yamaguchi model of the NN 
interaction, and Pnn stands for the PEST model of the NN interaction. The first column indicates 
the reference on which the resonant i]N model is based and the second column indicates the r/N 
scattering length (in fm) of that model. 



Ref. 




YtinYnn 


Y v nPnn 


RqNYNN 


RrjNPNN 


[6] 


0.72+10.26 


2.38+13.04 


2.59+12.37 


2.43+11.60 


2.46+H.62 


[7] 


0.75+10.27 


2.38+13.41 


2.70+12.64 


2.56+11.65 


2.61+H.72 


[8](D) 


0.83+i0.27 


2.57+i4.51 


3.21+i3.34 


3.00+U.87 


3.10+i2.03 


[8] (A) 


0.87+10.27 


2.54+i5.20 


3.47+13.78 


3.20+11.94 


3.36+12.19 


[8](B) 


1.05+10.27 


0.06+18.66 


4.27+17.05 


4.20+i2.24 


4.81+13.19 


[8](C) 


1.07+i0.26 


-0.42+i9.21 


4.53+i7.60 


4.21+i2.04 


5.02+i3.14 
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